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ABSTRACT 

In an interference network, joint power and admission control aims 
to support a maximum number of links at their specified signal to 
interference plus noise ratio (SINR) targets while using a minimum 
total transmission power. In our previous work, we formulated the 
joint control problem as a sparse £q -minimization problem and re- 
laxed it to a i\ -minimization problem. In this work, we propose to 
approximate the £n -optimization problem by a p norm minimization 
problem where < p < 1, since intuitively p norm will approximate 
norm better than 1 norm. We first show that the £ p -minimization 
problem is strongly NP-hard and then derive a reformulation of it 
such that the well developed interior-point algorithms can be applied 
to solve it. The solution to the £ v -minimization problem can effi- 
ciently guide the link's removals (deflation). Numerical simulations 
show the proposed heuristic outperforms the existing algorithms. 

1. INTRODUCTION 

Power control is an effective tool for interference management in 
cellular, ad-hoc, and cognitive underlay networks fl"ll2ll3ll4ll5l |6ll71 
8 |. The prevailing formulation of power control aims to use a mini- 
mum total transmission power to support all links in an interference 
network at their desired SINR targets. A longstanding issue associ- 
ated with power control is that the problem often becomes infeasible, 
i.e., it is not possible to simultaneously support all links in the net- 
work at their SINR targets. In this case, we must adopt a joint power 
and admission control approach to selectively remove some links 
from the network so that the remaining ones can be simultaneously 
supported at their desired SINR levels. Our goal is to maximize the 
number of simultaneously supportable links at their required SINR 
targets while using a minimum total transmission power. 

Theoretically, the joint power and admission control problem is 
known to be NP-hard to solve to global optimality 1 1 3 1 and to ap- 
proximate to constant ratio global optimality |6|, so various heuris- 
tic algorithms [12 3 , 6 , 7 . 8 1 have been proposed for this problem. 
Among them, the reference 1 1 1 proposed a convex approximation- 
based algorithm for the joint power and admission control problem. 
Instead of directly solving the original NP-hard problem, the basic 
idea of the proposed linear programming deflation (LPD) algorithm 
in (T) is to approximate the problem by an appropriate convex prob- 
lem. The solution to the approximation problem can be used to check 



the feasibility of the original problem and guide link's removals. The 
removal procedure is terminated until all the remaining links in the 
network are simultaneously supportable. The recent work 16] devel- 
oped another LP approximation-based new linear programming de- 
flation (NLPD) algorithm for the joint power and admission control 
problem. In (6), the joint power and admission control problem is 
first equivalently reformulated as a sparse to -minimization problem 
and then its £i -convex approximation is used to derive a LP, which 
is different from the one in (T]. Again, the solution to the derived 
LP can guide an iterative link removal procedure, and the removal 
procedure is terminated if all the remaining links in the network are 
simultaneously supportable. 

Based on the sparse £q -minimization reformulation in |6|, this 
paper proposes a new deflation algorithm based on p (0 < p < 1) 
norm minimization for the joint power and admission control prob- 
lem. Compared to the £\ -minimization problem, the p norm mini- 
mization problem is closer to the original £q -optimization problem. 
The ^-approximation problem is solved by applying the efficient 
interior-point algorithm in 1 9 1 to solve its equivalent reformulation. 
Numerical results show that the proposed algorithm compares favor- 
ably with the existing approaches 1 1 , 2 , 6 1 in terms of the number of 
supported links, the total transmission power, and the CPU time. 

Notations: We adopt the following notations in this paper. We 
denote the index set {1, 2, ■ ■ • , K} by M Lowercase boldface and 
uppercase boldface are used for vectors and matrices, respectively. 
For a given vector x, the notations max{x}, [x]fc and ||x|| p := 
^2 k |[x]fcj p (0 < p < 1 J] stand for its maximum entry, its fc-th 
entry, and its p norm, respectively. In particular, when p = 0, ||x||o 
stands for the number of nonzero entries in x. Finally, we use e to 
represent the vector of an appropriate size with all components be- 
ing one and I to represent the identity matrix of an appropriate size, 
respectively. 

2. PROBLEM FORMULATION 

Consider a A'-link (a link corresponds to a transmitter-receiver pair) 
interference channel with channel gains gkj > (from transmitter j 
to receiver fc), noise power r/k > 0, SINR target 7^ > 0, and power 

'Strictly speaking, ||x|| p with < p < 1 is not a norm, since it does 
not satisfy the triangle inequality. However, we still call it p norm for conve- 
nience in this paper. 



budget pk > for k, j 6 AC := {1, 2, • • • , if}. Denote the power 
allocation vector by p = (pi,P2, ■ ■ ■ ,Pk) T and the power budget 
vector by p = (pi, pa, ■ ■ ■ ,px) T - Treating interference as noise, 
we can write the SINR at the fc-th receiver as 



SINR*. = 



gkkPk 



Vfc £ /C. 



The joint power and admission control problem can be math- 
ematically formulated as a two-stage optimization problem [T]. 
Specifically, the first stage maximizes the number of admitted links: 



max |<S| 
p, s 



(1) 



s.t. SINR fe > 7 fc , k e S C /C, 
< p < p. 

We use So to denote the optimal solution for problem {T} and call it 
maximum admissible set. Notice that the solution for l(T]( might not 
be unique. The second stage minimizes the total transmission power 
required to support the admitted links in So: 



s.t. SINR fc > 7fc, fc G So, 
< Pfc < Pk, k £ So- 



(2) 



Due to the special choice of So, power control problem Q is feasi- 
ble and can be efficiently and distributively solved by the Foschini- 
Miljanic algorithm [4|. 



where < a < ai := l/e T p, and [q e ]fc measures the excess 
transmission power 1 1| that the transmitter of link k needs in the nor- 
malized channel in order to be served with its desired SINR target 
(assuming all other links keep their transmission powers unchanged). 
Notice that the formulation ([3j is capable of finding the maximum 
admissible set with minimum total transmission power. Since prob- 
lem {3} is still NP-hard, we further consider its i\ -convex approxi- 
mation (equivalent to a LP) 



min 

qe, q 

s.t. 



||q e ||i +ap q 

q e = c Aq, 
< q < e. 



(4) 



By solving {4|, we know whether all links in the network can be 
simultaneously supported or not. If not, we drop the link 



fc = argmax I ^ (ja fcj | [q^ + \a jk \ [q e ] fc J 

An easy-to-check necessary condition 

(/x + ) T e - (//" + e) T c > 



(5) 



(6) 



for all links in the network to be simultaneously supported is also 
derived in (6), where fi + — max{/i,0} , fj, — max{— fJ.,0} , 
and /j, = A T e. The necessary condition allows us to iteratively 
remove strong interfering links from the network. In particular, we 
remove the link ko according to the scheme 



3. REVIEW OF THE NLPD ALGORITHM 

Since the developed algorithm in this paper follows the similar idea 
as the NLPD algorithm in [6|, we first briefly review the NLPD al- 
gorithm in this section. The basic idea of the NLPD algorithm is 
to update the power and check whether all links can be supported 
or not. If the answer is yes, then terminate the algorithm; else drop 
one link from the network and update the power again. The above 
process is repeated until all the remaining links are supported. 

We begin with the introduction of an equivalent normalized 
channel on which the NLPD algorithm is based. In particular, we 
use q = <?2, • • ■ , Qk) T with qk = Pk/pk to denote the normal- 
ized power allocation vector, and use c = (ci, C2, • • ■ , ck) T with 
Cfc = (7fe?7fe) / (gkkpk) > to denote the normalized noise vector. 
We denote the normalized channel matrix by A £ R K x K with the 
(fc, i)-th entry 



akj 



IkQkjPj 

gkkPk 



if k = j; 

if Mi- 



In fact, \akj\ is the normalized channel gain. It is simple to check 
that SINR fc > 7 fc if and only if [Aq - c] fc > 0. 

In (6), we reformulate the two-stage joint power and admission 
control problem l|T]l and {2| as a single-stage optimization problem 



mm 
q E , q 

s.t. 



||q e ||o + ap q 

q e = c Aq, 
< q < e, 



(3) 



k = arg max < ^ | a kj | + ^ | a jk | + 



Cfc 



(V) 



until |6} becomes true. 

The NLPD algorithm can be described as follows. 



The NLPD Algorithm 
Step 1. Initialization: Input data (A, c, p) . 

Step 2. Preprocessing: Remove link ko iteratively according to 
{7| until condition {6| holds true. 

Step 3. Power control: Solve problem (H); check whether all links 
are supported: if yes, go to Step 5; else go to Step 4. 

Step 4. Admission control: Remove link fco according to set 
K = K./ {k } , and go to Step 3. 

Step 5. Postprocessing: Check the removed links for possible 
admission. 

4. A P NORM MINIMIZATION DEFLATION ALGORITHM 

In this section, we develop a new deflation algorithm based on £ p - 
minimization for the joint control problem {TJ and l(2j. As seen in 
Section [5] the original £o -minimization problem is successively 
approximated by the i\ -minimization problem |4j in the NLPD al- 
gorithm. Intuitively, the p (0 < p < 1) norm minimization problem 



mm 
q e . q 

s.t. 



||q e || p + ap T q 

q c = c - Aq, 
< q < e 



(8) 



should approximate {3} better than j4). This is the motivation for 
the development of the new deflation algorithm based on the p norm 
minimization for the joint power and admission control problem. 

Comparing the 1 norm minimization problem |4| and the p norm 
minimization problem ([8j, we see that problem Q is convex while 
problem |8]( is nonconvex (for its objective function is nonconvex 
with respect to q e ). Generally speaking, convex problems are rela- 
tively easy to solve, while nonconvex optimization problems are dif- 
ficult to solve. However, not all nonconvex problems are hard since 
the lack of convexity may be due to an inappropriate formulation. 
In fact, many nonconvex optimization problems admit a convex re- 
formulation. Therefore, convexity is useful but unreliable to test the 
computational intractability of an optimization problem. A more ro- 
bust tool is the computational complexity theory 1 10 1. 

We show that problem ([8} is strongly NP-hard. The proof is 
based on a polynomial time reduction from the MAX-2UNANIMITY 
problem, which is shown to be strongly NP-hard in II II . 

Theorem 4.1 The £ p -minimization problem {8} is strongly NP-hard 
ifO<P< 1- 



The complexity result in Theorem |4. 1 1 motivates us to approx- 
imately solve problem {8}. Next, we first give a reformulation of 
problem (8}, and then propose to use the interior-point algorithm de- 
veloped in 1 9 1 to solve it. 

Theorem 4.2 The l p -minimization problem ([8} can be equivalently 
reformulated as 



min E fc [q e ]fe + ap T q 
s.t. q e = c - Aq, 

< q < e, q e > 0. 



(9) 



A rigorous proof of Theorem |4.2| shall be given in the journal ver- 
sion. Here we just shed some light on why it does not harm optimal- 
ity to restrict q e > 0. Notice that problem |9j is equivalent to 



min 

qe, q 

s.t. 



||q e || p + ap T q 

q e = c Aq, 

< q < e, q e > 0. 



Thus, to show the equivalence of <j8j and (9}, it suffices to show 
that any optimal solution (q e , q) of ([8} always satisfies q e = c 
Aq > 0. In fact, assume the contrary that \ fZ + \ > 1, where K. + = 
{ k | [q e ]fe > }, IC = — { k j [q e ]fc = }. Then by the Balancing 
Lemma (see [6 Lemma 1]), we can appropriately reduce the power 
of links in JC + U so that both the first term and the second term 
in the objective of ([8} are strictly decreased. 

Based on Theorem |4. 2 1 by introducing a slack variable s > to 
problem |9), we see that problem ([S| is actually equivalent to 



_min_ J2 k Ml + a P T <l 



is, q, f 
s.t. 



q e = c — Aq, s + q = e 
q > 0, q e > 0, s > 0. 



(10) 



Now, we can apply the interior-point algorithm in |9| to solve prob- 
lem (To). Similar to (9], we can prove that the potential reduction 



interior-point algorithm returns an e-KKT 1 12] or e-global solution 
of problem (IP) in no more than O f ( — — ) log(-)^ iterations. 

One may ask why we wish to use interior-point algorithms to 
solve problem {TO}? The reasons are the following. First, the objec- 
tive function of problem (10) is differentiable in the interior feasible 
region. Moreover, we are actually interested in finding a sparse solu- 
tion q e of problem <jT0j>; if we start from a solution, some of whose 
entries are already zero, then it is very hard to make it nonzero. In 
contrast, if we start from an interior point, the interior-point algo- 
rithm may generate a sequence of interior points that bypasses solu- 
tions with the wrong zero supporting set and converges to the true 
one. This is exactly the idea of the interior-point algorithm devel- 
oped in 1 12 1 for the nonconvex quadratic programming. 

The proposed p norm minimization deflation (PNMD) algorithm 
is given as follows. There are two unclear points in the PNMD al- 
gorithm. One is how to compute the parameter a in problem (lOj, 
and the other is which removal strategy will be used in the admission 
control step. Next, we make clear of these two points, i.e., we shall 
use a given in (12) and the removal strategy (13) in the new deflation 
algorithm. 

In the NLPD algorithm |6|, the parameter a is given by 



ciai, if p(I — A) > 1, 

c 2 min{ai, a 2 } , if p(I — A) < 1, 



(11) 



where < c\ < C2 < 1 are two constants, and ot\ is determined by 
the equivalence between problem {3} and the joint problem {I} and 
(|2j, and Q2 is determined by the so-called "Never-Over-Removal" 
property. Since the £ p -minimization problem l[8j is closer to the to- 
minimization problem (3), we relax the parameter a in \ \\ \ to 



ciQi, if p(I — A) > 1, 

min {c2«i, 0302} , if p(I — A) < 1, 



(12) 



where C3 > C2, < ci, C2 < 1 are three constants. 

Having obtained the solution (q e , q, s) of problem < [ 1 0| >, we use 
the removal strategy called SMART rule in (3 1 to drop the link ko 
according to 



feo = argmax < ^ \a-kj\qj + ^ \ajk\qk + c k 



(13) 



The above operation can be interpreted as removing the link with the 
largest interference plus noise footprint in the normalized network. 



The PNMD Algorithm 
Step 1. Initialization: Input data (A, c, p) . 
Step 2. Preprocessing: Remove link ko iteratively according to 
(7} until condition |6| holds true. 

Step 3. Power control: Compute the parameter a and solve prob- 
lem (10) ; check whether all links are supported: if yes, go to Step 
5; else go to Step 4. 

Step 4. Admission control: Remove link ko according to some 

removal strategy, set K, = K,/ {&o} , and g° to Step 3. 

Step 5. Postprocessing: Check the removed links for possible ad- 
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Fig. 1. Average number of supported links versus the number of 
total links. 
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Fig. 2. Average CPU time versus the number of total links. 

5. NUMERICAL SIMULATIONS 

We generate the same channel parameters as in (T| in our numeri- 
cal simulations, i.e., each transmitter's location obeys the uniform 
distribution over a 2 Km x 2 Km square and the location of its cor- 
responding receiver is uniformly generated in a disc with radius 400 
m; channel gains are given by gkj = l/^fcj (VJfe, j £ /C), where dkj 
is the Euclidean distance from the link of transmitter j to the link of 
receiver k. Each link's SINR target is set to be 7^ = 2 dB (V k £ K.) 
and the noise power is set to be r)k = —90 dBm (V A; £ K,). The 
power budget of the link of transmitter k is pk = 2p™ m (V k £ AC), 
where p™ 111 is the minimum power needed for link k to meet its SINR 
requirement in the absence of any interference from other links. 

The parameter p in problem ([HJ is set to be 0.5 and the ones 
in ( |12[ > are set to be c\ = C2 = 0.2 and C3 = 4. The number of 
supported links, the total transmission power, and the CPU time are 
the metrics we employ to compare the performance of the proposed 
PNMD algorithm with that of the LPD algorithm in Q, the Algo- 
rithm II-B in |2], and the NLPD algorithm in (6|. All figures are 
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Fig. 3. Average transmission power versus the number of total links. 

obtained by averaging over 200 Monte-Carlo runs. 

Figs.[T]|2] and[3]indicate that the PNMD algorithm can take less 
CPU time to support more links while with less total transmission 
power than the existing algorithms (except the Algorithm II-B). As 
shown in Fig. [3] the Algorithm II-B transmits the least power among 
the tested algorithms. This is because the Algorithm II-B supports 
the least number of links; see Fig. [T] In particular, compared to 
the NLPD algorithnj^] the proposed PNMD algorithm can support 
(slightly) more links with much less total transmission power, and at 
the same time takes less CPU time. 

The performance improvement of the proposed PNMD algo- 
rithm over the NLPD algorithm is mainly attributed to the l v - 
approximation problem (JSJ. The simulation results in Fig. [T]and 
Fig. [3] show that the admissible set S\ obtained by the proposed 
PNMD algorithm based on the l p -approximation problem ([8J is 
"better" than the admissible set 52 obtained by the NLPD algorithm 
based on the £\ -approximation problem l|4j, i.e., although the cardi- 
nality of the two admissible sets 5i and 52 is nearly equal to each 
other, it takes much less total transmission power to support the 
links in 5i than to support the links in 52. This is consistent with 
our intuition that the p (0 < p < 1) norm minimization problem (8J 
is capable of approximating the £0 -minimization problem {3} better 
than the £\ -minimization problem |4} and the fact that the maximum 
admissible set for the joint power and admission control problem 
may not be unique. 

6. CONCLUSIONS 

In this paper, we have developed a p (0 < p < 1) norm minimiza- 
tion deflation algorithm for the joint power and admission control 
problem. Numerical simulations show the proposed algorithm out- 
performs state-of-the-arts in [T][2][6) m terms of the number of sup- 
ported links, the total transmission power, and the CPU time . 

2 To the best of our knowledge, the NLPD algorithm is so far the best 
removal-based algorithm for the joint power and admission control problem. 
It is shown in [6] that the NLPD algorithm can achieve more than 98% of 
global optimality in terms of the number of supported links when K < 18. 
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